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VCG  MLASUREMENT  AND  DISPLAY 


INTRODUCTION 


Aircrewmen  referred  to  the  USAF  School  of  Aerospace  Medicine  have 
a vectorcardiogram  (VCG)  taken  as  part  of  their  routine  evaluation.  The 
X,  y,  and  z lead  signals  are  recorded  as  frequency-modulated  (FM)  signals 
on  analog  magnetic  tape,  and  conventional  Polaroid  photographs  are  taken 
of  the  VCG  loops.  These  data  are  then  made  available  to  research  cardi- 
ologists through  a system  which  digitizes,  filters,  and  measures  the 
data  to  provide  parameters  that  previously  were  manually  derived.  The 
digitizing  and  filtering  portion  of  the  system  ha.s  been  documented  (11). 

This  report  defines  the  computer  algorithms  that  provide  VCG 
measurements  for  the  cardiologist.  These  algorithms  require  VCGs  t'nat 
have  been  baseline  adjusted  and  converted  to  millivolts.  Most  of  these 
measurements  are  the  ones  normally  derived  by  cardiologists  for  clinical 
evaluation.  Some  additional  measurements,  which  are  more  difficult  for 
the  human  to  obtain,  are  provided  to  assist  in  the  cardiologist's 
evaluation  and  for  investigational  purposes;  e.g.,  identifying  specific 
features  to  use  in  an  automated  VCG  analysis. 

The  conventional  Polaroid  photographs  provided  by  thf  VCG  equipment 
used  at  USAFSAM  (6)  display  data  at  2.5-msec  intervals,  or  400  samples/ 
sec.  The  data  acquired  by  the  digitizing  system  are  sampled  every 
2 msec,  or  500  samples/sec.  To  provide  displays  comparable  to  the  pho- 
tographs, the  digitized  data  are  linearly  interpolated  to  provide  an 
effective  sampling  rate  of  400  samples/sec.  In  the  Polaroid  photographs, 
each  sample  is  presented  as  a teardrop  to  assist  in  identification  and 
to  determine  direction  of  motion.  Our  system  represents  each  inter- 
polated sample  as  a plus  sign  which  increases  in  size  with  time. 


VECTOR  DEFINITIONS 


The  computer  algorithms  use  several  vectors  and  their  length.  For 
clarity  and  accuracy,  the  following  convention  will  be  used.  When 
vectors  are  referenced,  the  symbols  will  be  under'*  ined;  when  the  vector 
length  is  referenced,  the  symbol  will  not  be  underlined.  For  example, 
the  heart  vector,  h(t) , whose  components  are  the  x,  y,  and  z lead  signals 
for  time  (t): 


h(t) 


'x(t) 

y(t) 

z(t) 


and  Its  length,  h(t): 
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h(t) 


it) 


h(t) 


(t)‘ 


y(t)" 


z(t)' 


The  time  derivative  of  the  heart  vector  (the  velocity  vector)  and  its 
length  are  defined  as  VV(t)  and  VV(t),  where 


VV(t) 


d h(t)/dt  = 


’dx(t ) /dt 
dy (t)/dt 
dz(t)/dt 


and 


VV(t) 


=v/y 


VV  (t)  VV(t) 


= \/(dx(t 


)/dt)^  + (dy(t)/dt)"  + (dz(t)/dt)‘ 


Figure  I shows  the  velocity  vector,  VV(t),  partitioned  into  two 
components:  one,  VVj^j^(t),  parallel  to  the  heart  vector;  and  the  other, 

VVj(t),  perpendicular  to  it.  These  two  components  represent  the  de- 
rivative of  the  heart  vector  resulting  from  stretching  and  swinging. 

When  the  heart  vector  remains  constant  in  length  and  rotates, 

yVj^^(t)  = 0.  When  the  heart  vector  remains  fixed  ivi  space  and  lengthens 

or  shortens,  VVj (t)  = 0.  Finally,  it  can  be  shown  that  the  length  of 
VVj^^(t),  which  we  designate  W^j^(t),  is  the  time  derivative  of  the 
length  of  the  heart  vector: 


W^^(t)  = dh(t)/dt 
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The  Incremental  heart  vector 
from  one  sample  to  the  next. 


is  the  change  in  the  heart  vector 


Sh(k)  = h(k)  - h(k-l) 


x(k)  - x(k-l) 
y(k)  - y(k-i) 
z(k)  - z(k-l) 


where  k signifies  the  kth  sample.  Ali  velocity  vectors  are  approxi- 
mated by  incremental  vectors,  wiiich  means  that  the  velocity  vector 
approximations  are  measured  in  mV/sample  (2.5  msec)  rather  than  mV/msec. 


Two  mean  QRS  vectors  are  measured.  The  first  is  the  usual  (stan- 
dard) mean  QRS  vector,  QRS  , whose  components  are  obtained  from  the 

mean 

individual  signals  (4).  In  the  x and  y leads,  the  usual  mean  QRS  values 
are  defined  as  the  algebraic  sum  of  the  amplitudes  of  the  signal  at  the 
minimum  of  the  Q wave  and  maximum  of  the  R wave  (see  Fig.  2).  In  the 
(-z)  lead,  the  usual  mean  QRS  values  are  defined  as  the  sum  of  the 
amplitudes  of  the  signal  at  the  maximum  of  the  R wave  and  the  minimum 
of  the  S wave  (see  Fig.  3).  If  any  of  the  waves  do  not  exist  in  a par- 
ticular lead,  the  value  of  the  amplitude  is  assumed  to  be  zero.  After 
the  components  are  determined,  the  usual  mean  QRS  vector  is  assembled. 


QRS 

mean 


QRS 

(x) 

mean 

QRS 

(y) 

mean 

QRS 

(2) 

mean 

J 

The  second  mean  QRS  vector  measured  is  designated  the  true  mean 
QRS  vector  (1,  4).  It  is  a vector  along  the  vector  mean  of  all  the 
heart  vectors  comprising  the  QRS  loop. 


QRS 


mean 


x(k) 

y(k) 

2(k) 


where  N and  N are  the  indexes  at  which  the  QRS  loop  begins  and  ends, 

Q J 

respectively.  Figure  4 shows  an  example  of  the  true  mean  QRS  vector. 

The  true  mean  QRS  vector  can  also  be  viewed  in  another  manner.  Con- 
sider the  X lead  shown  in  Figure  5.  The  area  under  the  QRS  wave  can 
be  approximated  by  the  area  enclosed  in  the  rectangles  (7). 
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Figure  2. 


Definitions  used  for  the  usual 
mean  QRS  values  in  x and  y leads 
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Figure  3.  Definition  used  for  the  usual 
mean  QRS  values  in  the  z lead. 
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Tills  can  Ih> 


N 


N 


A = 


j:  x(k) 

k=N.. 


AT  = 


AT  I x(k) 
k=N^ 


where  AT  is  the  width  of  the  rectangles  (sampling  interval).  Observe 
that  rectangles  below  the  zero  line  have  negative  areas.  The  average 
value  of  X during  the  ORS  loop  can  be  approximated  by 


AT 


i:  x(k)/(N 


Nq  + 1 ) AT 


Treating  leads  y and  z in  the  same  manner  and  assembling  the  results 
into  a vector  yields  the  equation  for  the  true  mean  QRS  vector.  There- 
fore the  true  mean  vector  is  based  on  area,  while  the  usual  mean  vector 
is  based  merely  on  the  extreme  amplitudes  of  the  Q,  R,  and  S waves. 


The  vectors  measured  are  presented  to  the  cardiologists  as  projec- 
tions into  four  planes — frontal,  left-sagittal,  transverse,  and  eigen. 
When  the  values  of  these  planar  projections  are  identified  for  examina- 
tion. polar  coordinates  will  be  used.  The  references  for  and  positive 
directions  of  the  angles  are  given  in  Figure  6.  The  angles  are  always 
less  than  180°  (+  270°  is  presented  as  its  equivalent  - 90°,  etc.) 


x(lef t) 
0° 


I’RONTAl, 

f 90*' 
y (inferior) 


Figure  6.  Examples 


0 z (posterior)  z (posterior) 
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of  angular  measurement  in  the  various  planes. 
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POINT  DEFINITIONS  AND  EXAMPLES 


The  first  task  of  the  measurement  portion  of  the  routine  is  to 
identify  one  heart  beat.  This  is  performed  by  locating  a fiducial 
point  on  each  R wave.  This  fiducial  point  is  designated  as  the  Npth 
sample  and  is  used  to  define  a window  of  Interest  1 sec  or  the  mean 
RR  interval,  whichever  is  shorter,  in  length.  This  window  is  posi- 
tioned so  that  45%  of  it  occurs  before  Np  and  55%  occurs  after.  These 
percentages  are  empirically  chosen  after  examination  of  more  than 
1,000  vectorcardiograms.  This  fiducial  point  is  defined  as  the  sample 
corresponding  to  the  maximum  velocity  vector. 

The  X,  y,  and  z lead  data  have  been  previously  processed  to  define 
a zero  reference  and  adjust  all  data  to  it  (12);  however,  the  zero- 
reference  definition  in  the  preprocessor  is  a weighted  sum  of  the  TP  and 
PR  segments.  A decision  to  use  the  PR  segment  as  the  zero-reference  was 
made  because  USAFSAM  cardiologists  plan  to  use  these  algorithms  with 
data  obtained  at  higher  heart  rates,  such  as  those  during  exercise. 
Therefore,  this  zero  reference  is  shifted  from  the  commonly  used  TP 
segment.  Since  the  PR  segment  is  sloped  in  many  cases,  a more  precise 
definition  is  required  for  the  amplitude  measurement  reference  point. 
After  examining  more  than  1,000  VCGs,  a decision  was  made  to  define  the 
PR  segment  as  follows.  The  150-msec  interval  (61  samples)  of  the  veloc- 
ity vector  data  prior  to  the  fiducial  is  searched  to  determine  the  15- 
msec  segment  of  least  activity.  The  average  value  of  the  heart  vector 
for  that  15-msec  segment  is  defined  as  the  zero-reference  value.  This 
region  of  least  activity  is  defined  as  the  15-msec  segment  of  data  for 
which  the  average  velocity  vector  is  a minimum  (see  Fig.  7). 


Figure  7.  Determination  of  the  15-rasec  baseline  interval  (Np^) . 
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The  computer  algoritlim  for  determining  tlie  15-msec  zero-reference 
segment  finds  a minimum  for 
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To  agree  with  the  measurements  obtained  by  the  cardiologists,  the  origin 
of  the  Frank  coordinate  system  is  adjusted  to  this  zero  reference  by 


x(k)‘ 

x(k)  •'  Xj^j 

y(k) 

= 

y(k)  - y^^^ 

z(k) 

z(k)  - 

where  N is  the  number  of  samples  in  the  beat. 

D 

Similar  logic  is  used  to  determine  the  beginning  and  ending  points 
of  the  three  loops  (P,  QRS,  and  T)  examined  by  cardiologists.  First  a 
region  of  interest  is  defined  for  each: 

(1)  For  the  P loop,  the  time  interval  (N__)  from  the  beginning  of 
the  beat  to  the  beginning  of  the  PR  segment.  ^ 

(2)  For  the  QRS  loop,  the  time  Interval  from  the  end  of  the  PR 
segment  to  100  msec  after  the  fiducial  point. 

(3)  Only  the  end  of  the  T loop  is  determined  (its  beginning  is 
defined  to  be  the  end  of  the  QRS  loop)  from  the  time  interval  between 
the  end  of  the  QRS  loop  anc  the  end  of  the  beat. 

In  each  region  of  interest,  the  maximum  velocity-vector  length,  VV, 
is  determined,  and  a threshold  is  defined  as  a predetermined  percent 
of  this  maximum.  A 10-msec  window  (5  samples)  is  moved  forward  from 
the  maximum  W until  all  points  in  it  are  less  than  the  threshold.  A 
second  10-msec  window  is  moved  backward  in  time  until  all  VV  values  are 
less  than  the  threshold  (see  Fig.  8).  The  beginning  and  ending  of  the 
loops  are  defined  as  the  rightmost  point  in  the  backward-moving  window 
and  the  leftmost  in  the  forward-moving  window. 
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klL 


Figure  8.  An  example  of  the  loop-finding  algorithm. 

This  loop-defining  algorithm  works  well  in  most  cases;  however, 
there  are  exceptions.  To  handle  some  of  these,  the  algorithm  includes 
the  following  additional  rules. 


P Loop 

If  the  duration  of  the  P wave  is  greater  than  150  msec,  the  thresh 
old  is  raised  by  10%  of  the  maximum  VV  and  the  process  is  repeated.  If 
necessary,  this  incrementing  is  continued  until  the  criteria  are  satis- 
fied or  the  threshold  reaches  50%  of  the  maximum  VV.  If  the  latter 
occurs,  a flag  is  set  and  processing  continues  for  the  next  measurement 


QRS  Loop 


If  the  threshold  criteria  are  not  satisfied  for  the  beginning  of 
the  QRS  loop,  the  beginning  is  defined  as  the  end  of  the  zero-reference 
segment  (Np^  + 6th  sample).  Likewise,  if  the  criteria  for  the  ending 

are  not  satisfied  within  100  msec  after  the  fiducial  point,  the  ending 

is  defined  as  100  msec  after  the  fiducial  (N  + 40th  sample). 

r 


T Loop 

If  the  threshold  criteria  are  not  satisfied  for  the  ending  window, 
the  threshold  is  raised  in  increments  of  10%  of  the  maximum  W,  and  the 
procedure  is  repeated  until  either  the  criteria  are  satisfied  or  the 
threshold  reaches  50%  of  the  maximum  W.  In  the  latter  case,  a flag  is 
set  and  processing  continues. 
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All  graphic  displays  are  generated  on  a Calcomp  Piotter.  Figure 

y is  the  first  display,  consisting  of  the  tiirec  original  x,  y,  and  z 

leads  plus  VV.  The  indexes  corresponding  to  the  beginning  and  ending 

of  tlie  P loop  are  defined  as  N and  N ; the  zero-reference  segment, 

PB  PE 

Np^;  and  the  beginning  and  ending  of  the  QRS  loop,  and  (for  J- 
junction).  Finally,  the  T loop  begins  at  the  J-junction  (Nj)  and  ends 
at  the  N th  sample. 

Next,  investigators  at  USAFSAM  (1,  7,  9)  are  interested  in  examin- 
ing the  VCG  in  its  ovna  plane;  therefore,  the  VCG  is  transformed  into  a 
coordinate  system  where  the  maximum  energy  is  contained  in  two  of  the 
signals  (the  eigenpiane) . The  basis  vectors  for  this  coorainate  system 
are  the  eigenvectors  of  the  matrix  (2,  3,  5,  10,  11,  13,  14). 


Z 

i=l 


'x.  y , z , ' 
11  1 
1 i 


Because  the  eigenvectors'  directions  are  not  unique,  and  in  order  to 
agree  with  the  work  of  the  others  (7-9,  11),  the  basis  vectors  are 
rotated  about  the  out-of-eigenplane  vector  so  that  the  final  basis 
vectors  (U , V,  and  W)  have  the  following  orientations.  _V  lies  along 
the  projection  of  the  true  mean  QRS  vector  (half  area)  in  the  eigen- 
plane,  and  W is  perpendicular  to  the  eigenpiane,  oriented  so  that  the 
rotation  of  the  QRS  loop  in  the  eigenpiane  (about  W)  is  positive  as 
defined  by  the  right-hand  rule.  Finally  W is  chosen  so  that  IJ,  V^,  and 
W are  pointed  upward  and  form  the  polar  vector  (7-9).  Also,  the  premise 
is  that  the  V lead  signal  optimizes  repolarization  evaluation  since  it 
lies  along  the  electrical  axis  of  tlie  heart.  The  individual  eigenlead 
signals  are  presented  to  the  cardiologist  as  in  Figure  10.  The  temporal 
points,  previously  determined,  are  also  indicated  on  these  plots.  The 
()RS  and  T loops  in  the  eigenpiane  are  displayed  as  shown  in  Figure  11. 


The  traditional  planar  tracings  are  displayed  to  the  card iolog.ist  as 
shown  in  Figures  12,  13,  ard  14.  In  these  figures,  the  frontal  plane  is 
viewing  the  patient  from  the  front  (in  the  z direction),  the  sagittal 
plane  is  viewing  the  patiert  from  his  left  side  (in  the  minus  x direction), 
tlie  transverse  plane  is  viewing  the  patient  from  the  head  downward  (in  the 
y direction).  Kach  of  these  planar  projections  is  scaled  so  that  the  loop 
being  displayed  fills  a 5-i  m square.  The  scale  factor  is  a reference 
ail  justed  as  illustr.ited  in  Figures  11-14.  This  is  indicated  by  the  size 
ot  the  "+"  symbols;  as  time  increases,  the  symbols  become  larger.  The 
interval  between  these  symbols  is  2.5  msec.  Also  shown  on  these  planar 
views  i.s  the  location  (indicated  by  an  asterisk)  of  the  vectors  prior  to 
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Figure  12.  Planar  tracings  of  frontal,  left  sagittal, 
and  transverse  QRS  loops. 
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Figure  14.  Planar  tracings  of  frontal,  left 
sagittal,  and  transverse  P loops. 
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Finally,  three  vectors  of  importance  to  USAFSAM  cardiologists  are 
displayed  as  shown  in  Figure  15.  These  are  (1)  a vector  designated  J — 
the  J-junctlon,  h(Nj);  (2)  a vector  designated  M — the  heart  vector  40 

msec  (16  samples)  after  the  J-junction,  li(Nj+16);  and  (3)  a vector 

designated  T,  which  is  0.1  mV  in  length  and  with  the  same  orientation 
as  the  maximum  T vector.  From  the  first  two  vectors  (J^  and  M) , both 
ST  depression  and  ST  slope  can  be  determined.  The  ^ unit  vector  is 
plotted  to  provide  an  indication  of  the  relative  positions  of  M and  ^ 
to  the  maximum  T. 
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Figure  15.  Planar  J,  M,  and  T vectors. 


Printed  Output 

Along  with  the  graphic  displays  described  in  the  previous  sections, 
three  pages  of  computer-printed  results  are  presented  to  the  cardiologist 
Figures  16  and  17  are  examples  of  the  first  two  pages — the  basic  data. 

On  the  first  line,  the  VCG  number  identifies  the  patient.  Although  the 
computer  algorithms  were  designed  to  be  used  on  average  beats  (12),  they 
can  also  be  used  on  relatively  nolsefree  individual  beats.  In  each  case. 
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the  particular  beat  examined — 1st,  2il,  3d,  or  average — Is  identified 

on  the  first  line  of  the  printout;  next,  the  patient's  heart  rate  is  given, 
based  on  the  RR  interval.  The  second  line  contains  P,  QRS,  T,  and  QT 
duration  times.  These  are  derived  from  the  temporal  points  previously 
de  f i ned . 


QT  duration 
T duration 
QRS  duration 
P duration 


(2 . 5 msec) 
(2.5  msec) 
(2.5  msec) 
(2.5  msec) 


"‘tf.  - V 
‘“te  - “j' 


(Nj  - y 


The  percent  of  out-of-eigenplane  energy  (3d-line  display)  is 
determined  by 

X . / (X  + X . , + X . ) 
min  max  mid  min 


where  the  X’s  are  the  eigenvalues  used  in  determining  the  eigenplane, 
and  the  subscripts  indicate  their  relative  size.  On  the  4th  line,  the 
matrix  for  transforming  the  vectors  from  the  Prank  coordinate  system  to 
the  eigenplane  coordinate  system  is  printed  as  a unit  vector  and  an 
angle  (theta).  The  unit  vector  is  the  vector  about  which  the  Frank 
coordinate  system  can  be  rotated  into  alignment  with  the  eigenplane 
system,  and  theta  is  the  required  rotation  angle. 


A series  of  measurements  is  made  on  the  preprocessed  x,  y,  z,  u, 
and  V leads.  The  R-wave  and  T-wave  maximum  amplitudes  and  their 
ratios  are  given  for  each  lead.  The  components  of  the  heart  vector 
at  the  J-junction  along  with  the  components  of  the  slope  of  the  ST 
vector  are  also  measured. 

ST  , = [h  (N,  + 16)  - h (NJ]  /40 

— slope  — J — J 

The  units  of  the  slope  are  mV/sec.  The  remaining  measurements  pertain 
to  the  vector  loops. 


A series  of  measurements  is  made  that  allow  the  cardiologist  to 
examine  the  terminal  forces  (vectors)  of  the  QRS  loop.  These  are  the 
heart  vectors  at  5,  10,  20,  30,  and  40  msec  into  the  loop,  and  30,  20, 
and  10  msec  from  the  end  of  the  loop.  To  quantify  how  fast  the  vectors 
are  moving,  the  average  value  of  the  lengths  of  the  projections  of  the 
incremental  vector  into  the  four  planes  is  determined  for  (1)  the  initial 

30  msec  of  the  QRS  loop,  (2)  the  time  segment  between  initial  and  final 

30  msec  of  the  QRS  loop,  and  (3)  the  final  30  nsec  of  the  QRS  loop. 

These  measurements— labeled  MEAN  ACT.  1ST,  MEAN  ACT.  MID,  AND  MEAN  ACT. 

LAST — are  estimates  of  the  mean  magnitude  of  the  projections  of  the 
velocity  vector  onto  the  four  planes,  where  units  are  mV/sample. 
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DIBECTIOM 


Figure  17.  Page  2 of  printout. 


Measurements  labeled  YQ-ZR  WIDTH  and  YQ-ZR  HEIGHT  relate  to  the 
QKS  loop  and  are  defined  as  follows.  In  the  frontal  plane,  YQ  width  is 
the  maximal-initial-left-superior  force,  the  value  of  x at  the  end  of 
the  Q wave  in  y (y  = 0);  in  the  left-sagittal  plane,  the  maximal-initial- 
superior-posterior  force,  the  value  of  z at  the  end  of  the  Q wave  in 
y (y  = 0);  in  the  transverse  plane,  the  maximal-initial-lef t-poster ior 
lorce,  the  value  of  z at  the  end  of  Q wave  in  x (x  = 0) ; in  the  eigen- 
plane,  the  value  of  u at  the  end  of  the  Q wave  in  v (v  = 0) . In  the 
frontal  plane,  height  is  the  maximal-superior-initial  force  (YQ  height); 
in  the  eigenplane,  the  maximum  amplitude  of  the  Q wave  in  v.  Examples  of 
these  measurements  are  shown  in  Figure  18. 
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Figure  18.  Examples  of  the  height  and  width  measurements 
in  various  planes. 


21 


The  J,  M,  and  T vector  lengtlis  and  orientations  are  presented 
(Figs.  16,  17)  with  the  following  exception.  Instead  of  giving  the 
angular  location  of  the  .1  and  M veciors  in  the  planes,  only  the  quad- 
rant in  which  they  are  located  is  given.  For  example,  RS  signified 
rig.htwanl  and  suiierior.  Also,  ani’les  hetwei-n  the  projections  of  ^ 

.and  ^ onto  the  four  planes  are  given.  The  reference  for  raeasuring  the 
angles  is  the  ^ projection,  and  positive  rotations  are  defined  in 
Figure  6. 

The  usual  (STD)  and  maximum  QRS  vectors'  lengths  and  orientations 
are  printed  (Fig.  16)  with  only  the  orientation  of  the  true  mean  QRS 
vector  given. 

The  values  of  the  projections  of  the  heart  vector  at  the  peak  of 
QRS,  T,  and  P loops  in  the  four  planes  are  given,  where  the  peaks  of 
the  loops  are  defined  for  each  plane  as  the  point  at  winch  the  projec- 
tion of  the  heart  vector  onto  that  plane  has  its  maximum  length.  These 
projection  vectors  are  labeled  MAX. QRS,  MAX.T,  and  MAX.P,  respectively. 
Also  given  are  the  ratio  of  the  length  of  MAX.T  to  MAX. QRS,  along  with 
the  angle  between  these  vectors.  The  reference  for  this  angle  is  the 
MiXX.QRS  vector,  and  tlie  definition  of  positive  angle  is  given  in 
Figure  6. 

The  last  (3d)  sheet  of  the  printout  (Fig.  19)  is  used  to  evaluate 
the  program  and  to  record  the  cardiologist's  diagnosis.  The  evaluation 
section  is  self-explanatory.  USAFSAM  cardiologists  enter  codes  for  the 
following;  scalar  codes  (A),  vector  codes  (B),  SSI  code  (C — locally 
used  diagnostic  classification),  and  the  T.C.D.  (D — terminal  conduction 
delay  in  msec).  These  are  used  with  the  performance  evaluation  of  an 
automated  VCG  analysis  program. 


APPLICATIONS 

The  previous  sections  have  described  the  logic  of  VCG  measurement 
algorithms  that  are  being  applied  to  a set  of  vectorcardiograms  used  in 
developing  an  automated  VCG  analysis  program  specifically  designed  for 
the  USAF  Central  ECG  Repository.  Host  of  these  measurement  algorithms 
are  applicable  for  routine  clinical  use  in  any  automated  VCG  evaluation 
program.  Several  of  these  algorithms  (e.g.,  the  determination  of  the 
eigenplane  and  its  spatial  orientation)  were  developed  for  investiga- 
tions designed  to  Identify  potentially  unique  diagnostic  features  of 
the  VCG.  The  algorithms  can  be  readily  modified  for  widespread  applica- 
tion in  many  electrocardiographic  analysis  programs. 
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Figure  19.  Cardiologist  review  sheet  (page  3 of  printout). 
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